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Abstract 

We prove nonlinear stability for finite amplitude perturbations of plane Couette flow. 
A bound of the solution of the resolvent equation in the unstable complex half-plane 
is used to estimate the solution of the full nonlinear problem. The result is a lower 
bound, including Reynolds number dependence, of the threshold amplitude below 
which all perturbations are stable. Our result is an improvement of the corresponding 
bound derived in P). 

1 Introduction 

Plane Couette flow is stable to infinitesimal perturbations for all Reynolds numbers 
Finite amplitude perturbations on the other hand can induce transition to turbulence. 
Thus there is a threshold amplitude below which all perturbations decay eventually (there 
may be transient growth) and above which some perturbations lead to turbulent flow. 
In the present paper we derive a bound depending on the Reynolds number R, for this 
threshold. 

The stability of shear flows, and plane Couette flow in particular, has been extensively 
studied, we refer to El and ^! for more background information. Early work focused 
on the eigenvalues of the linearized problem. The question of the asymptotic value of the 
threshold amplitude for large R was starting to be investigated in the early 1990s 
A dependence 0{R~'^) with 7 k. 5/4 is supported by computations in [Sj. The asymptotic 
analysis of gives the result 7 ~ 1. Another approach to this question is found in [H] 
where an estimate of the resolvent equation is used to prove an upper bound for the 
exponent, 7 < 21/4. In this paper we use the same approach and sharpen the bound of 7 
obtained in This is possible because we work with a norm which weighs the different 
coefficients of the velocity vector, the weights depend on the Reynolds number. It is 
numerically demonstrated in jSj that, in the weighted norm, the resolvent grows linearly 
with i?, as opposed to quadratically which is the case of the energy norm used in [3]. 
Because of the weighting the exponent 7 will be different for the different components of 
the velocity vector. The precise result is stated in Theorem ^ 
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The proof of the result of this paper is carried out with the same techniques as in P]. 
However, the starting point is a different resolvent estimate. Also, the energy norm fits 
the structure of the Navier-Stokes equations, the modified norm we use does not. This 
leads to several technical complications. 

In Section 121 we give the mathematical formulation of the problem, the Navier-Stokes 
equations, the linearized equation and the resolvent equation. Then we state the main 
result in Theorem ^ The proof of the theorem is organized as follows. In Section we 
analyze the linear problem. Using the result of mentioned above, we derive an estimate 
for higher derivatives of the solution of the resolvent equation (this solution is the Laplace 
transform of the velocity field). By inverse transformation we obtain a time dependent 
estimate. In Section ^ the linear estimate is used to derive the threshold value for the 
nonlinear problem. 

In Appendix^ we have collected the definition of the weighted norm and other non- 
standard norms used in the paper. Technical results including estimates of norms of the 
nonlinearity are collected in Appendix El In Appendix we prove a theorem which will 
be applied to show that the perturbation tends to zero if the suitable a priori estimate 
holds. The result of AppendixlHl explains the choice of left hand side in the linear estimate 
of Theorem |21 in Section |21 



2 The main result 

We choose the coordinate system so that the velocity field of Couette flow is given by 



in the domain 

= {a; G : -1< X2 < 1} . 

We use bold letters to denote vectors and subscripts to identify the components. The 
functions and pco = const constitute a stationary solution of the nondimensionalized 
Navier-Stokes equations 

ut + G{u) + gradp = -^Au, 

divii = 0. (2.1) 
Here R denotes the Reynolds number and G is the following nonlinear differential operator. 

k=l 

We have no-slip boundary conditions, it = on dQ, and the initial condition 
u{x,0) = uc^{x) +v^^\x), 
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where 

^(0) is 

the initial perturbation. We will assume that v^^^ E H*. Local existence of 
a classical solution of H2.1() is proven in |5j (theorem 5 together with theorem 7 on p. 161 
and p. 167 respectively of [5:). Furthermore we will derive a priori estimates which (for 
sufficiently small initial data) allows extension of the local solution to a global solution by 
successive application of the above result. 

We investigate the stability of Couette flow by linearization of the Navier-Stokes equa- 
tions. We denote the perturbation by v which thus is related to u according to 

U = Uco+ V. 

Below we will apply the Laplace transform and for this purpose we want a problem with 
homogeneous initial data. This is accomplished with the introduction of w by 

= e-*i;(°) + w. (2.2) 

The function w satisfies the following problem 

wt + gradp = Lw + e"*(L + l)v^^^ - G {e^^v^^^ + w 

diYW = 0, 

t(;(a;,0) = 0, (2.3) 

where 

L--A- — - 

R ^"^ dxi 

For w we also have no-slip boundary conditions. If we assume that and w are small 
and neglect G in (|2.3j) then we see that w satisfies a linear problem with forcing depending 
on the initial perturbation. This is the linearized equation for Couette flow. 

We will first consider the forcing in (|2.3j) as a given function, denoted by / = f{x,t) 
and derive estimates for the solution. Later we will return to the particular form of the 
forcing in (|2.3() and its dependence on the initial data and w. 

Application of the Laplace transform now give the resolvent equation 

sw + gradp = Lw + /, 
divw = 0, 

w = 0, xedn. (2.4) 

This problem was investigated in where the following estimate of the solution was 
obtained 

||w||^ < Ci?^||/|l^, Res>0. (2.5) 

See Appendix ^ equation (|A.1|) for the definition of the modified (m-)norm. We have 
collected the definitions of all norms used in this paper in Appendix Throughout 
the paper we will use C to denote constants which appear in inequalities, we do not use 
subscripts to identify different constants. We emphasize that C denotes an "absolute 
constant" and does not depend on the forcing, the initial data or the Reynolds number, 
which is the only parameter in the problem. 

We now state our main result in the following theorem. 
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Theorem 1. The resolvent estimate h2. implies that there is a constant 6 > such that 



if 



lim 

t— >oo 



\vk{-,o)\\Hi ^ = 1,3, 



Ik2(-, 0)11^4 < 

Comparing this theorem to the previously obtained bound on the threshold in we 
see that the i?-exponent for the first and third i7-components is improved from 5.25 (in [S]) 
to 3 and the exponent for the second component is improved from 5.25 to 4. 

It may appear excessive to require that the initial data is as smooth as H'^. The result 
in |n] require the iJ^-norm of the initial data to be small. However, in |^ a completely 
different approach is used with the application of semi-group methods and we are not 
aware of any results on the threshold amplitude obtained with this approach. It may 
be possible to weaken our smoothness requirement using an appropriate local existence 
theorem for the Navier-Stokes equations which incorporates the smoothing property. We 
also note the error in 0[ p. 193] where it is stated that small H'^-novm implies stability, 
for their result to hold small H^-norm is required. 

If one is interested in a particular perturbation v^^\x) = a(p{x), and the threshold 
coefficient a{R), then it is insignificant which HP-norm is required since 

3 The linear estimate 

In this section we derive estimates for the solution of the following problem in terms of 
the forcing function / 

Wt + gradp = Lw + /, 
divw = 0, 

w{x,0) = 0. (3.1) 

The results are collected in the inequality of the theorem below. In the statement of 
the theorem the Hi- and //2-iiorms appear for the first time, their definition is given in 
Appendix 1X1 

Theorem 2. The resolvent estimate \2. 5\) implies that the solution of equation \3. 1\) 
satisfies 



\W{;T)\\1^+I^ 





< C 



\w{;t)\\% + \\wt{;t)\\ljdt 

T 



+ CR' [ {\\f(.;t)fm + \\ft{;t)\t)dt. (3.2) 
Jo 
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For the remainder of the paper we wih suppress the notation "(•, t)" which indicates at 
which time a norm is evaluated, this will be clear from the context. 

In Section r3.1l we analyze the Laplace transformed problem (|2.4|) and then in SectionEIS 
we use these results to conclude the proof of Theorem |2 



3.1 Estimates of the transformed functions 

The starting point for this section is inequality (|2.5|1 which we will use to derive the 
following lemma. 



Lemma 1. If Res > 0, the solution of \2.4\ l satisfies 

\\w\\\<CRWm- 

The estimate is uniform in s. 

Proof. By partial integration, the following identity is easily derived 



(3.3) 



Re {w, Lw) 



1 ^ 

-y 

k=l 



dw 



+ Re {wi,W2). 



(3.4) 



Now, we take the L^-inner product of w and the resolvent equation (|2.4|) . Using (|3.4j) and 
the fact that w is solenoidal we obtain 



1 



Res||-u;f + - V 



k=l 



dw 


2 


dxk 


= Re 



{wi,W2) + {W, f) 



(3.5) 



The right hand side in ()3.5|) can be estimated using the inequality ab < + 6^/2e 

which holds for all e > 0. Without explicit mention we will use this inequality several 
times below. Here we take a = \\w\\, b = \\ f\\ and e = 1/R. We also use the fact that the 
L2-norm is smaller than the modified norm {R > 1), and obtain 



Re 



iwuw2) + {w,f) <CR\\f\\ 



I m' 



Combining (|3.6|) with (|3.5|) and using Re s > we have derived 



E 

k=i 



dw 



dxk 



< CR^Wfl 



(3.6) 



(3.7) 



It remains to estimate J2{'w) (defined in Appendix 0. To do this we differentiate the 
resolvent equation with respect to xi and X3. We give the details in the xi-case. We 
obtain 



dw dp dw df 

hgrad-— = L- h ^— . 

dxi dxi dxi dxi 



(3.8) 



Note that L commutes with d/dxi and d/dx^ but not d/dx2- We take the inner product 
of (|3.8p with dw/dxi and obtain 



Re s 



dw 



dxi 



k=l 



W 



dxidxh 



Re 



dwi dw2 



dxi ' dxi 



+ 



dw df 

dxi ' dxi 



(3.9) 
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The right hand side is estimated in the following way. 

' dwi dw2 ' 



Re 



dxi ' dxi 



+ 



< 



1 

4R 



dxf 



dw df \ 

dxi ' dxi J 

+ R\\w2f + 



1 



dxf 



+ R\\ff- 



(3.10) 



Now we insert the estimate (|3.1U() into ()3.9p and cancel the second derivatives in the right 
hand side with the corresponding terms on the left hand side. Further simplification leads 
to 



1 ^ 

-y 



k=l 



dxidxf^ 



<C{^\\w\t + R\\ffm] <CR\\f\ 



2 

m' 



The same procedure applies to the resolvent equation differentiated with respect to X3. 
We have proven 



J2{w)<CRWm- 

A combination of 1)2. 5() . 1)3. 7|) and 1)3. 11() gives the lemma. 
Remark. To obtain (|3.4j) we used 
{w,gvadp) = —{divw,p) = 0, 



(3.11) 



(3.12) 



where the first equality follows by partial integration and the second since w is solenoidal 
(divib = 0). It would be preferable to work with modified (m-)norms throughout the 
paper, the reason one cannot do this is 1)3. 12() . If one takes the modified inner product of w 
and the resolvent equation then it is not possible to eliminate the pressure using (|3.12|) . 
The modified inner product is defined by 

{u,v)m = iui,Vi) + R^{U2,V2) + (^Sj^^s)- 

This is the explanation of the simplifications with the use of the energy (L2) norm men- 
tioned in the introduction. 



3.2 Estimates of the time dependent functions 

Before we prove Lemma |21 below and Theorem [51 we state three inequalities which will be 
of use. The following is derived from equation (|3.1|) using (|3.4() 



d ,, ,,9 
dt" " R 



1 ^ 

-y 



k=l 



dw 



dxk 



< \iw,f)\ + \{wi,W2)\. 



(3.13) 



In the same way we obtain the next inequality, for ^ = 1, 3, starting from (|3.H) differentiated 
with respect to xi and X3 respectively 



d_ 

di 



dw 



dxi 



1 ^ 



fc=i 



w 



dxidx}- 



< 



dxf 



,f 



+ 



dxf 



W2 



(3.14) 
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We apply Plancherel's formula for the Laplace transform to (|3.3|) , with the imaginary axis 
as integration contour. This yields 



\w\\% dt < CR^ 

"Hi 



I fir dt 



(3.15) 



which holds for all T > 0, see 13| p. 235-239] for this application of Plancherel's formula. 
Lemma 2. The solution of IIS. 1\) satisfies 



dw 


2 


dw 


2 . 


dxi 


+ 


dx3 


< CR / 


Jo 



IfWmdt, 



(3.16) 



for all T > 0. 

Proof. From ()3.13|) we obtain 



\w\\^ < -WwWi + —11/ 



P +-\\wP 

^ II llm' 



e > 0. 



We integrate this from t = to T, choose e = 1/R and use 1)3. 15(1 . 



\wr < CR 



T 



\f\\jyi_dt. 



We estimate the right hand side in ()3.14() according to 



52 



dx 



Wi 

— ,W2 



+ 



d w 

dxf 



< 



dxf 



+ R\\w2\r + 



4R 



dxf 



+ R\\f\ 



dw 




<CR [ 


dxi 


Jo 



recall that / = 1,3. This inequality is inserted into 1)3. 14() . We move the second derivative 
terms to the left hand side and integrate in time to obtain 



This completes the proof of Lemma I 

Now we turn to the proof of Theorem [2 The function Wt satisfies the same PDE as w 
(with p and / replaced with pt and respectively) and inhomogeneous initial data, 

wt{x,0) = f{x,0) =: f^^\x). 

We use the symbol := to indicate that the expression on the "colon side" is defined in 
terms of the expression on the other side. Now we introduce u by 



wt = e-'f'^°^ +u. 

The function u satisfies the following problem. 

ut + gradpt = Lu + ft + e'^L + 1)/^, 
diva; = 0, 
u{x,id) = 0. 



(3.17) 
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The only difference between this problem and equation for w is the forcing. The in- 
equalities ()3.15() and l|3.16() can thus be applied here yielding estimates for u. Using (|3.17j) 
and the triangle inequality these can be converted to the following estimates for w. 



T 



IwtWi dt<C 

"Hi 



|/W|||^+i?l(L + l)/ 



(0)1 



tWm 



dt 



(3.18) 



and 




+ ^11(^ + 1)/ 



(0)1 



+ R 



turn 



dt 



/ = 1,3. (3.19) 



With ()3.15() and 1)3. 18(1 we have estimated the time integral in the left hand side of the 
inequality of Theorem El The i^2-norm at time T remains. To estimate the first order 
derivatives in the i72-norm we start from the following estimate which is easily obtained 
from (Tmil 



E 

k=l 



dw 



< CRi\iw,f)\ + \{w,wt)\ + \iwi,W2)\) . 



(3.20) 



Using the Cauchy-Schwartz inequality and (|3.16() and ()3.19|) we see that the right hand 
side of 1)3. 2U() can be estimated in terms of the right hand side of 1)3. 2|) . 
To estimate J2{w) we rearrange 1)3.14(1 and obtain 



k=l 



dxidxh 



< CR 



dwi dw-2 



dxi dxi 



+ 



dx^ 



,f 



+ 



dw dwt 

dxi ' dxi 



for I = 1,3. For the second term in the right hand side we have 



dxf 



< 



1 

4R 



dxf 



+ R\\f\ 



(3.21) 



(3.22) 



We insert ()3.22|) into ()3.21|) and move the second derivative term to the left hand side. 



The remaining two terms in the right hand side of ([3.21(1 are estimated using the Cauchy- 
Schwartz inequality and ()3.16|1 and (|3.19(1 . We see that J2{'w) can be estimated by the 
expression on the right hand side of the inequality of Theorem |21 This concludes the proof 
of the theorem. 



4 Derivation of the threshold bound 



Now we return to the full nonlinear problem (|2.3|1 for w. This is the same as ()3.1|1 if we 

take 




(0) 



' dxk ^ " dxk^ 



dxk 



+ Wk 



dw 

dxk 
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For simplicity we first consider 

f = e-HL + l)v^^^ + G{w). (4.1) 

The quadratic terms in v^^^ and the coupling terms between v^^^^ and w only add technical 
difficulties, the asymptotic value of the threshold amplitude is the same. To make the 
formulas below more lucid we introduce the following notation for expressions in the 
inequality of Theorem El 



N-iw,T) = \\w\\l^ + 1^ 



\w\\l^ + \\wt\\l^ 



dt, 



M{f,T) = CR^ 



\ff + 



l/llm "I" ll/tllm) 



Moif, T) =M{f,T)+C /(°) III, +RmL + l)/(0) II 



As long as the solution w exists, the estimate of Theorem |21 is valid with forcing ()4.1() . 
Expressed with our new notation we have 

N'iw,T) < Mo (^e-\L + 1)^(°) , r) + MiG{w),T). (4.2) 

Here we have used that w{x,0) = 0. The first term on the right hand side of 1)4.2(1 can be 
estimated in terms of v^^\ independently of T, as we show in AppendixEl inequality ()B.1|) . 
We have 



Mo (e-\L + l)v^^\T) <CR^\\v 



,(0)||2 



M. 



(4.3) 



Now we will determine the threshold value by assuming that w does not tend to zero. As 
a consequence of this assumption we will derive the inequality (|4.13|) (equivalent to ()4.14|) 
which bounds a norm of i?^'^^ from below. The inverse of this inequality gives the threshold 
which we thus prove by contradiction. 

After the above outline, we now turn to the details of the proof. We thus assume that 
||i(;(-, T)||oo does not tend to zero as T ^ oo. According to Theorem |S1 in Appendix IHl 
(with / = ||iu||j5^^) we then have 

lim J\f{w,T) = oo. 

T — ^oo 

In particular, there is a Tq such that 

M{w,To) = 2M (4.4) 

and 

J\f{w,T)<2M, T<To. (4.5) 

Using 1(4. 5() and the Sobolev inequality (|A.2|) we can now estimate the three terms of 
M{G{w),To) in terms of M. For the first term we have the inequality ()B.2|) from Ap- 
pendix IbI 



CR^\\G(w)f < CR^\ 



E 

k=l 



dw 

dxk 



(4.6) 
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The two factors in the right hand side can be estimate as follows 



|2 



\w\r^ < Cllwll^^ < CM, 



dw 



dxi 



< C\\w 



\% <CM. 

-Hi 



(4.7) 
(4.8) 



Inserting (|4.7|) and 1)4. 8() into 1)4. 6() . we obtain 
CR^\\G{w)\\^ < CR^'M^. 



(4.9) 



We use similar techniques and inequality HB.3|) and (|B.4|) for the remaining two terms, 



To 



\G{w)\\l^dt < CR^\\w\ 



To 3 

E 

k=l 



dw 



dxk 



dt < CR^M 



(4.10) 



and 



CR^ [ ° \\Gtiw)\\'idt < CR^V] 

^0 k=l 



dw 


2 j-To 


dw 


dxk 


m JO 


dt 



dt 



+ CR^\\w\ 



To 3 

E 

k=l 



dwt 



dxi 



dt < CR^M . 



(4.11) 



Combining (jCT7|) and (|OT|) we obtain 

M{G{w),To) < CR^JA^. (4.12) 
Now we insert (|0|) . (|^ and (|I?T^ into to obtain 

2M<M + CRV^^. (4.13) 
We divide the inequality by M., insert the expression 1)4. 3|) for M. and rearrange to 



,(0)| 



> 



1 



CR^' 



(4.14) 



We have thus shown that an initial perturbation which leads to instability must satis- 
fy ()4.14() . If the initial perturbation satisfies the inverse inequality then the instability 
assumption leads to a contradiction. In other words, if the initial perturbation satisfies 
the inverse inequality of 1)4. 14() then ||i<;(-, T)||oo as T tends to infinity. Noting the w 
is related to v by (|2.2|) . the proof of Theorem ^ is complete. 



A Definition and notation for norms 



In this paper we track the i?-dependence of all estimates we derive. This forces us to 
introduce some non-standard norms which weigh the different components of the vector 
function with coefficients depending on R. Also we need a Sobolev inequality (|A.2|) with 
a right hand side not containing all second derivatives. 
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We start, however, to introduce notation for the standard norms. The L^-norm for 
functions u : — > C is defined by 



u 



\u(x)\'^dx. 



n 



We use the same notation for the L^-norm of vector functions 
3 



u 



k=l 



It is clear from the argument if it is the scalar or vector norm which is intended. For the 
standard Sobolev spaces with "integration exponent" two we use the following notation 

\a\<k 



U 



where a = (ai, 02, as) is a multi- index. 

We refer to the weighted norms also as modified norms and use the subscript m to 
identify them. Let X denote any of the norms introduced so far, then the corresponding 
modified version is 



u 



\ui\\x + ^^||u2||l: + WusWx- 



In particular we have 



u 



\uif + R^\\u2f + \\u3f. 



(A.l) 

Before we give the Sobolev inequality HA.2|1 we define the semi-norm J2 consisting of 
selected second order derivatives. The choice of second order derivatives is dictated by the 
difficulty in estimating normal derivatives (derivatives with respect to the X2-coordinate) 



J2N 



u 



+ 



dxidx2 



+ 



dx2dx^ 



+ 



92 



u 



9x| 



Now we can define two similar norms 



u 



Hi 



U 



fc=i 



du 



and 



u 



H2 



k=l 



dxk 



du 



dxi 



+ J2{u) 



+ J2{u). 



The i/i-norm is greater than the //2-iiorm (i? > 1) for which we have the Sobolev inequali- 

ty 



\u\ 



<C\\u\\l^, 



(A.2) 



see jU p. 385]. By the left hand side in (|A.2|) we of course mean 

||u||oo = max ess sup |n^.(a;)[. 
fce{i,2,3} ccen 
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B Auxiliary estimates 

To bound M.o{e~^{L + 1)v^^\T), we must estimate the following four terms 



I = Ci?2||(L + l)w(0)||2, 

rT 



CB? f e-'%L + l)v('^\ldt = CR'\\{L + l)v('^\l=:II, 
Jo 



77I=||(L + iyo)|||^ 



and 

jy = II (L + 1)2^(0) II 

For the first term we have 



For the second term 



dxi 



H2- 



+ CR^v'f^? <CRVA 



if 2, 



The third and fourth terms are bounded according to 
III<CR^v^^'^\?„^ , IV KCRV^ml,^ . 

— II llii*,m' — II llii*,m 

Combining the above results we obtain 

Mo {e-\L + l)v^'\T)<CR'\\v^'^\l,^^. 
To bound M.{G{w),T) we need the following estimates of the nonlinearity 



\G(w) 



k=l 



fe=l 



|2 
loo 



dw 



dxk 



k=l 
2 



dw 



Wk- 



< 



k=l 
3 

E 

k=l 



dwi 



dxi 



dxk 

3 
k=l 



dw 



+ R^ 



dxk 



dxi 



+ 



Wk 



dxk 



\Wk\\lo 



dw 



dxi 



<\MIY1 



k=l 



dw 



dxk 



(B.l) 



(B.2) 



(B.3) 
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and 



\Gtiw) 



\m 



d I V - dw 



< 



E 

k=l 



dwk 


2 


dw 


dt 


oo 


dxk 



< 



dw 



dt 



2 3 



• k=i 



dw 



dxt 



k=l 
3 

+ E 

k=l 



dwk dw 



dt dxi 



+ 



dwt 



Wk 



dxk 



\Wk\\lo 



dwi 



dxk 



.E 

k=l 



dwi 



dxi 



(B.4) 



C A Sobolev inequality for asymptotic stability 



Let / : [0, oo) 



i. The theorem below gives a condition which imply that f{t) as 



the time {t) tends to infinity. In this paper we will choose f{t) = t) H^:^^ , combined 



with the inequality ()A.2|) . the theorem can then be used to prove that 
Theorem 3. // 



{\f{t)?+\f'it)\')dt< 



oo, 



0. 



(C.l) 



then 



lim f{t) = 0. 
Proof. The required Sobolev inequality is 

poo 

SUp|/(t)|2<C / (|/(t)|2 + |/'(t)|2)dt, 
t>T JT 

where C is independent of T. For a proof, see [31 appendix 3]. 
Clearly 

|/(r)|2<sup|/(t)p. 

t>T 

Because of ()C.1|) we also have 

lim / (l/(t)|2 + |/'(t)|2)dt = 0. 



(C.2) 



(C.3) 



Combining (|(y.2|) . ()('.3|) and ()('.4|) completes the proof. 



(C.4) 
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